Many-body correlations of electrostatically trapped dipolar excitons 
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We study the photoluminescence (PL) of a two-dimensional liquid of oriented dipolar excitons in 
Ina;Gai_a;As coupled double quantum wells confined to a microtrap. Generating excitons outside 
the trap and transferring them at lattice temperatures down to T = 240 mK into the trap we 
create cold quasi-equilibrium bosonic ensembles of some 1000 excitons with thermal de Broglie 
wavelengths exceeding the excitonic separation. With decreasing temperature and increasing density 
n < 5 X 10^"^^ we find an increasingly asymmetric PL lineshape with a sharpening blue edge and 
a broad red tail which we interpret to reflect correlated behavior mediated by dipolar interactions. 
From the PL intensity I{E) below the PL maximum at Eo we extract at T < 5K a distinct power 
law I{E) ~ (Eo — _E)^'"' with -\a\ ^ -0.8 in the range Eo — E oi 1.5-4 meV, comparable to the 
dipolar interaction energy. 



Weakly interacting bosons confined in an external po- 
tential and cooled to very low temperatures such that the 
thermal de Broglie wavelength becomes comparable to 
the inter-particle distance can form a Bose-Einstein con- 
densate (BEC). In this new state of matter, a large frac- 
tion of the bosons condense into the energetically lowest 
quantum state of the external potential, and form a cor- 
related state with a macroscopic wavefunction. BEC has 
been observed in different systems such as ultra cold di- 
luted atomic gases |TJ [5] or superfluids |3] . More recently, 
condensation of non-equilibrium quasi-particles in solids, 
namely cavity exciton polaritons [IHB] and magnons 7j, 
have been reported as other examples of BEC. In an effort 
to realize BEC of excitons in quasi-equilibrium, predicted 
already in the 1960s [H H] , coupled double quantum wells 
containing two-dimensional systems of long-living and 
spatially indirect excitons (IX) have been intensively ex- 
plored for more than a decade [see, e. g. P 1^ [rUHTBl [TSl - 
[20]], but fully convincing signatures of a BEC ground 
state are still missing. Such IX consist of an electron 
and a hole, spatially separated in the adjacent wells of 
coupled double quantum well (CDQW) heterostructures 
and bound by their Coulomb attraction. With their ori- 
ented dipolar nature they form a complex and rather 
strongly interacting bosonic systems with internal struc- 
ture and not yet fully understood ground state proper- 
ties. Only recently the influence of dipolar interactions 
on the ground state properties of such excitonic ensem- 
bles has been more intensely theoretically investigated 

To achieve a quantum phase transition, an efficient 
trapping of suitable IX densities thermalized to cryogenic 



temperatures is needed |15[ I23j . Going beyond previous 
experiments, we made special efforts to avoid a thermal 
imbalance between the temperature of the trapped IX 
and the lattice. In choosing the InGaAs material systems 
and resonantly exciting direct excitons in the CDQW at 
an energy below the band gap of the GaAs substrate we 
assure that most absorption of laser light occurs only in 
CDQW under the tight focus of the exciting laser. In ad- 
dition, our special micron scale trap design sketched in 
Fig. [ija) serves to generate IX on the slide gate several 
micrometer away from the trap and collect only precooled 
and neutral bosonic IX under the trap gate. The elec- 
trostatic trap uses the in-plane variation of the quantum 
confined Stark effect (QCSE) to spatially separate the 
generation of IX from the IX recombination in the trap. 
This allows to cool the device in a '^He refrigerator down 
to lattice temperatures of 240 mK and to study the pho- 
toluminescence (PL) of typically several 1000 trapped IX 
with high spatial, energetic, and temporal resolution. In 
the corresponding IX density regime and with decreasing 
temperature at which the thermal de Broglie wavelength 
becomes larger than the excitonic spacing, we observe 
the integrated intensity of the IX line to increase, while 
the linewidth narrows and a characteristic lineshape with 
a sharp blue edge develops. Since this happens at tem- 
peratures and IX densities, at which the kinetic energy 
of the dipolar excitons falls well below their repulsive in- 
terdipolar energies of typically 2 meV, we interpret the 
characteristic change in the PL lineshape and intensity 
as a signature of a correlated liquid of bosonic IX as the- 
oretically discussed in Ref. [S] . Increasing the IX density 
at lowest temperature by increasing the power of the in- 
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cident laser we also observe an initial decrease of the PL 
linewidth combined with an increase of the PL intensity 
and lineshape asymmetry. Analyzing the PL lineshape 
I{E) in the energy range Eq — E <1.5-4meV below the 
PL maximum at Eq we find a characteristic temperature- 
independent power law I{E) ^ (Eq — i?)^'"' , discernable 
at temperatures T < 5 K and moderate densities and 
indicating an edge-like singularity at Eq. 
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FIG. 1: (Color online) (a) Sketch of the trapping device 
and detection scheme, incorporating the projection of a scan- 
ning electron micrograph of the trap (T) (with diameter of 
6/im), slide (S), and two guard gates (G). Voltages Vr, Vs 
and Vg are applied with respect to the back contact (green), 
(b) Schematic CDQW band diagram. Ec and denote con- 
duction and valence band. Electron and hole wavefunction of 
ground states g'^°™"='^QW ^^uppcrQw indicated in green 
and red, respectively, (c) Pseudo-3D picture of excitonic PL 
from the trap with spatial and energetic resolution [Vg = 
0.70 V, Vt = 0.35 V, and Vs = 0.45 V, TLattico = 4K, PLasor 
— 8.0 fiW, i?LasGr — 1.494 eV]. (d) Corresponding ground state 
energies of electrons in the lower QW (g'^o^^'^Q^) a^d heavy 



holes in the upper QW (/i/i" 



uppGrQW\ 



versus y-position across 



the trap. The IX trapping potential l/jx is shown containing 
a hybridized dipolar liquid. 



In our device two 7nm thick Ino.iiGag.sgAs quantum 
wells (QW), separated by a 10 nm thick GaAs barrier, 
are embedded between a n-doped GaAs back contact and 
6 nm thick semi-transparent titanium gate electrodes de- 
fined by e-beam lithography [Fig. [IJa)] |2]. Voltages ap- 
plied to the gates define an electric field F perpendicular 
to the QW plane, and cause a red shift of the IX energy 
of AEix = — pF, with | p |= ed and d the electron-hole 



distance [Fig. [T|b)]. In-plane variation of F thus guides 
IX from the generation spot, where IX are generated via 
resonantly pumped direct intra well excitons, into the 
tunable trap |24j . Employing a confocal microscope with 
two objectives we can separate exciton generation and 
PL detection, measured in transmission, at different lo- 
cations [23]. An image of the PL intensity in Fig. [IJc) 
demonstrates the trapping efficiency. As sketched in Fig. 
[T|d), the trapping potential Vj^ is attractive for neutral 
dipolar IX but it repels unbound electrons. The low mo- 
bility of holes at low temperatures combined with them 
being trapped by the in-plane disorder potential let us as- 
sume that holes are essentially localized. In contrast, the 
wavefunctions of lighter electrons, confined in the lower 
QW and bound to holes localized in the upper QW as 
IX, might overlap and hybridize as their effective Bohr 
radius Vc ~ 20 nm is comparable to half the excitonic sep- 
aration dix (<iix/2 ~ 28 nm at nix = rio = 4 x 10^'^-^) 
[2¥ and the QW separation d « 17 nm, as indicated in 
Fig. [l|d). Caused by the repulsive dipolar interaction of 
the bound IX in the hybridized system one may expect 
correlated many-body behavior, similar to a Fermi edge 
singularity ^5^. 

The trap collects a cold IX ensemble of varying den- 
sity and a typical lifetime of about 200 ns caused by the 
reduced overlap of the electron and hole wavefunctions 
forming the IX. To explore the exciton thermodynamics, 
we analyze the temperature dependence of the IX PL 
collected from the center of the trap. With decreasing 
temperature an increasingly asymmetric PL lineshape is 
observed as exemplarily shown in Fig. [2f^a) . As displayed 
in Fig. ^h) , we also observe a strong increase of the inte- 
grated PL intensity accompanied by a narrowing of the 
PL linewidth. From these experimental observations, we 
deduce that many-body correlations start to become im- 
portant with decreasing temperature, finally resulting in 
an edge-like singularity. To achieve a deeper insight into 
the underlying exciton dynamics, we analyze the high 
(blue) and low (red) energy sides of the PL line sepa- 
rately. From the blue side of the PL lineshape we extract 
an apparent exciton gas temperature T^^pp inside the trap 
|24) . To this end, we assume only temperature-induced 
broadening and neglect all other broadening, e.g., by dis- 
order [Fig.[2jc)]. Using local resonant generation of direct 
excitons in the CDQW we minimize radiative heating 
|24) . Since the thermal equilibration via phonons hap- 
pens on time scales below 1 ns, much shorter than the IX 
lifetime, we estimate the real IX temperature in the trap 
to be much closer to TLattico than to Tlipp. The satura- 
tion of Tlipp ~ 3 K at low temperatures, corresponding to 
energy of ^ 0.3 meV, is likely to be caused by inhomoge- 
neous line broadening and temporal fluctuations. 

Many-body processes can be observed in characteris- 
tic PL lineshapes |25J HHj . From recently published the- 
oretical [57 and experimental works fT^ on lineshapes 
in quantum structures, one can conclude that different 
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many-body phenomena with distinct energy scales are re- 
flected by different power laws in the lineshape of discrete 
optical transitions at energy E. Such power laws corre- 
spond to different temporal dynamics, becoming faster 
with increasing Eq — E [57], where Eq denotes the en- 
ergy of the PL maximum. In previous PL experiments, 
the red tail of the spectra was attributed to the random 
disorder potential [Sll EQ] or the interplay of disorder po- 
tential and BEC [H]. As in Refs. [IlllTj, we analyze 
the red tail of the PL lineshape by plotting the PL in- 
tensity logarithmically versus the logarithm of the PL 
energy in a log(/) vs. log(i?o — E) diagram [Fig. [2]jd)]. 
In the log-log display it is possible to identify power laws 
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FIG. 2: (Color online) (a) Energy-resolved PL spectra, for 
different lattice temperatures, detected from the center of tlie 
trap [Vg = 0.70 V, Vt = 0.32 V, and Vs = 0.35 V, PLascr = 
2.7 fiW, -ELaser = 1.4289 cV]. (b) Corresponding integrated PL 
intensity (left axis) and full linewidth at half maximum (right 
axis), (c) Exciton gas temperature Tapp extracted from the 
blue tail of the PL line [PLascr = 1.4 /iW]. (d) Red tail of PL 
normalized logarithmically plotted as a function of Eo — E, 
with Eq the energy of the PL maximum [TLattico = 250 mK] . 
Color-coded (e) corresponding exponents and (f) intersection 
energies Eq — Ei at which the power laws in (d) change versus 
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FIG. 3: (Color online) (a) PL spectra, for different laser 
powers, detected from center of the trap. PL intensity is 
normalized to PLaser and the spectra are energetically shifted 
such that the Fermi-edge-like steepest slope in the PL vs. E 
spectra is at i?E ~ (original spectra see [21] )• (b) Integrated 
PL intensity plotted as a function of blue shift AS and corre- 
sponding density. Vertical lines in (b) and (c) are color-coded 
to reflect the different laser powers in (a), (c) Second expo- 
nent a (left axis) and linewidth (right axis) as a function of 
Prascr. (d) An upper limit is shown for Tkpp [Vg — 0.70 V, Vt 
= 0.32 V, Vs = 0.35 V, TLattico = 243 mK, ^Lascr = 1.4289 eV]. 



I{E) ~ {Eq— E)~^"-^ with different exponents a separated 
by kinks at which different power laws intersect p4j. In 
Fig. [2]Jd), and (e), we identify three different exponents 
associated with characteristic energy intervals [Fig. [2][f )] 
related to specific interactions [23] . The appearance of a 
few exponents in an optical correlation function seems to 
be typical for many-body correlated systems [27]. The 
most interesting exponent is the second exponent show- 
ing values of about -0.8. We relate this to a correlated 
low temperature many body phase of a two-dimensional 
dipolar liquid of parallel oriented excitons as theoretically 
predicted in Ref . [5^ . This power law rather abruptly dis- 
appears at TLattico ^ 5K [Fig. [2je)]. As predicted in Ref. 
[5] a correlated dipolar liquid is expected only at suf- 
ficiently low temperature and high enough trapped IX 
densities and should exhibits a rather gradual transition. 

We further investigate the PL spectra as a function of 
the exciton density, varied by changing the exciting laser 
power. In Fig. [sj^a), the PL intensity, normalized to the 
laser power, is plotted versus energy relative to the blue 
edge energy E-^ at which the PL intensity drops to half of 
the maximum value (original spectra see |24j ) . Around 
the PL maximum, the intensity rises more steeply than 
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laser power PLaser- Using the blue shift Ai? of the PL 
maximum with increasing laser power as a measure of the 
exciton density, we obtain the density dependence of the 
integrated PL intensity shown in Fig. [Sf^b) [53] . We point 
out that only a part of the nonlinear increasing PL inten- 
sity can be explained by a decreasing IX lifetime [2] ■ For 
-Ftascr < 1.2 /iW, we observe a decreasing linewidth with 
rising exciton density [Fig. [sjjc)]. A related observation 
in quantum wells was recently associated with screening 
of the disorder potential by repulsive IX-IX interaction 
pOj . The density dependent decrease of the linewidth 
is accompanied by finding the exponent -jaj « -0.8 in 
the log-log plot of the red tail [Fig. [sjjc)] . In the same 
density regime we find at low temperatures the exponent 
-|q| « -0.8 in the temperature-dependent spectra [Fig. 
[2]je)] which supports the our that this power law is asso- 
ciated with a correlated dipolar liquid. With increasing 
exciton density the dipolar interaction energies of the IX 
sum up to the observed blue shift AE [24] . The lineshape 
asymmetry thus can be interpreted as a broadened sum 
of individual excitonic transitions. In related studies on 
electrostatic traps fabricated on the same InGaAs het- 
erostructure with an effective trap diameter of down to 
about 100 nm and thus a trap area of roughly 1000 times 
smaller than those considered here, we observe discrete 
individual lines of linewidth below 0.5 meV which depend 
on laser power as characteristic for single excitons, biexci- 
tons and triexcitons, respectively. Further increasing the 
trap population, these merge at low temperatures into 
an asymmetric lineshape thus lending additional support 
to the interpretation given here in terms of correlated 
behavior. At high laser powers, the extracted exponent 
a starts to oscillate due to an interference substructure 
overlaying the PL lineshape [Fig.jsjjc), for details see [23]] 
and we see an expected homogeneous broadening |19j of 
the PL line caused by exciton-exciton scattering with in- 
creasing exciton density. From the blue side of the PL 
lineshape we extract a rising upper limit for the apparent 
exciton gas temperature Tapp [Fig. [sj^d)]. 

In conclusion, at temperatures and densities at which 
the thermal de Broglie wavelength (« 75 nm at 2K) is 
larger than the interexcitonic distance dix ~ 57 nm [24] 
at n = 4 X 10^° we observe a characteristic tem- 
perature and density-dependent transition to a specific 
quantum state of the IX ensemble. In this regime, dix 
is comparable to twice the effective excitonic Bohr ra- 
dius Te ~ 20 nm [24] as well as the QW separation {d w 
17nm). From the appearance of the strongly increasing 
integrated PL intensity as well as a PL linewidth decreas- 
ing with increasing density, we conclude that this state 
reflects a correlated many-body state mediated by dipole- 
dipole interactions. In the analysis of the red tail of the 
PL lineshape we observe at TLattice ^ 5 K a characteris- 
tic power law dependence of the PL intensity reflected in 
exponent « -0.8, associated with an energy between 
1.5 meV and 4 meV. This we interpret as a generic fin- 



gerprint of dipolar many-body correlations. A theoreti- 
cal understanding of this power law cannot be given at 
present but this very interesting observation in our opin- 
ion should motivate further theoretical and experimental 
studies in the field. 
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SUPPLEMENTAL MATERIAL: MANY-BODY 
CORRELATIONS OF ELECTROSTATICALLY 
TRAPPED INDIRECT EXCITONS 



FIG. SI: Cryogenic setup. Schematic diagram and photo- 
graph of our low temperature microscope with two diffraction 
limited confocal objectives embedded in a ^He refrigerator. 



CRYOGENIC SETUP 



RESONANT EXCITON EXCITATION 



In Fig. SI a scheme as well as a photograph of our 
experimental setup are displayed. The single mode fiber- 
based low-temperature microscope is composed of two 
confocal objectives. The microscope is embedded in a 
■^He refrigerator with a base temperature below 240 mK. 
The cryogenic setup exhibits an excellent long-time sta- 
bility in position against temperature sweeps and mag- 
netic field variations. Both objectives are diffraction lim- 
ited and thus have a focus spot size with a diameter 
comparable to the wavelength and can be individually 
positioned with a low-temperature piezo-ccyz-positioning 
unit ^4 . The positioning works at temperatures down to 
240 mK along all three space coordinates over millimeter 
distances and with a nanometer accuracy. The objective 
on top of the sample serves to generate excitons with laser 
radiation. The second objective below the sample is used 
to collect the photoluminescence (PL) light in transmis- 
sion and analyze it outside the cryostat, e.g. in a grating 
spectrometer. In order to efficiently cool the device chip, 
it is mounted on a cold gold plate and bonded several 
times to a cold ground. 



In our trapping configuration the exciton generation 
occurs on the slide gate typically spatially separated by 
9 /im from the center of the trap and the exciton trap 
is filled only with indirect excitons, pre-cooled to lattice 
temperature [5]. To further reduce the radiative heating 
in the sample and the excitation of free charge carriers, 
we resonantly pump the intra well direct exciton with a 
tunable diode laser. The photon energy of the excita- 
tion laser is below the GaAs band gap and, apart from 
some absorption in the semitransparent metal gates, can 
be only absorbed in the coupled double quantum wells. 
Accordingly, only the energy of the direct to indirect ex- 
citon conversion is deposited in the phonon bath of the 



sample. In Fig. S2 a photoluminescence excitation (PLE) 
measurement at TLattico = 246 mK is shown. Here, the 
integrated PL intensity of the indirect excitons is plotted 
as a function of the laser excitation energy exciting di- 
rect excitons on an unstructured gate 10 /im away from 
the detection spot. In Fig. [S2] we observe a strong IX 
PL when the direct exciton is resonant with the tunable 
laser energy and so direct excitons are efficiently pumped 
transforming into the detected indirect excitons. 
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the PL vs. E are E-e, « 0. The corresponding original 
spectra are shown in Fig. S3|b). To obtain these spec- 
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FIG. S2: Photoluminescence excitation measurement. 

Integrated PL intensity centered at 1.399 eV of the indirect 
excitons as a function of the laser excitation energy [l/batc = 

0.35V, PLaser(£LaBcr) = 4.75AtW, TLatticc = 246 mK] . 



DENSITY DEPENDENCE OF THE 
CORRELATED EXCITONIC 
PHOTOLUMINESCENCE 



To explore the tunability of the exciton density in the 
trap with illumination intensity, we vary the incident 
laser power focused on the slide quasi-statically by tun- 
ing the power transmitted into the microscope with an 
acousto-optical modulator. Thus it is possible to change 
the excitation power by three orders of magnitude. Fig. 
pSf^a) is identical to Fig. 3(a) in the main text and plots 
PL spectra, for different laser powers, detected from the 
center of the trap. The PL intensity is normalized to the 
excitation laser power and the spectra are energetically 
shifted such that the Fermi-edge-like steepest slopes in 



tra we integrate 55 seconds with a detection rate of 12.5 
counts per second in the PL maximum at 1.0 /iW. 



CONVERSION FROM BLUE SHIFT TO 
EXCITON DENSITY 



The quantum confined Stark effect causes a red shift 
of the excitonic energy of Ai^ix = — pF where p is the 
IX dipolc moment and F the electric field perpendicular 
to the QW plane. In a trapped IX ensemble all indi- 
vidual dipoles p = ed are aligned and perpendicular to 
the QW plane. The effective electric field, seen by the 
individual dipole moment of a single IX, is the sum of 
the externally applied electric field and the depolarizing 
dipole field of all other IX in the vicinity. With increas- 
ing exciton density the effective electric field decreases 
and causes a blue shift /S.E. For not too small densities 
n oi n ^ 2.5 x 10^° the simplest approximation is 
that the blue shift is proportional to the exciton density 
AE' = ^^^^n (see black fine in Fig. S4) where e is the 
dielectric constant. In the low-density limit the IX- IX 



interaction was theoretically investigated in Refs. [3H5] 
and results in a density and temperature-dependent cor- 
rection factor /(T, n) in Ai; = ^^nf{T, n). From this 
formula we find a nonlinear relation between blue shift 
and density as plotted by the red curve in Fig. |S4| The 
latter is used to convert blue shift AE to density in Fig. 
3(b) of the main text. 

A third naive T — approach to convert the mea- 
sured blue shift into a trapped IX density is to sum up 




E-E^ (meV) PL Energy (eV) pi_ Bi^e shift AE (meV) 



FIG. S3: Laser power dependence of the PL lineshape 
normalized versus original spectra, (a) Energy resolved 
PL spectra for different laser excitation powers (identical to 
Fig. 3(a) in the main text). The PL intensity is normalized 
to the excitation laser power and the spectra are energeti- 
cally shifted to correct the density caused blue shift. In (b) 
the corresponding original spectra are shown [Vg ~ 0.70 V, 
Vt = 0.32 V, and Vs = 0.35 V, TLatticc = 243 mK, ^Lascr = 
1.4289 eV]. 



FIG. S4: Conversion from blue shift to exciton density. 

The black straight line shows the linear conversion between 
the blue shift AE of the PL maximum energy and the cor- 
responding IX density. In the red curve a density dependent 
correction factor is included providing a better approxima- 
tion of the real exciton density. The upper green curve show 
the conversion between blue shift and trapped IX density for 
a two-dimensional hexagonal ordered crystal of parallel ori- 
ented dipolar excitons. 
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the dipole fields of the trapped indirect exciton ensemble 
acting on a single dipolar exciton. For this method it 
is assumed that the dipolar exciton liquid is ordered in 
a two-dimensional hexagonal crystal of dipoles arranged 
in parallel. A single dipolar exciton sees a blue shift, 
caused by another dipolar exciton in a distance rfix-ix, 

2 

ofW— —73^ where p — edis the IX dipole moment re- 

■^^ix-ix 

suiting from the electron-hole separation by a distance d 
in the two adjacent quantum wells and e is the dielectric 
constant. The total blue shift AE for a single exciton 
with a next neighbor distance dix is calculated by sum- 
ming up all contributions W of every trapped exciton. 
The next neighbor distance djx is converted by n = 
in the corresponding trapped indirect exciton density n. 



The green curve in Fig. S4 shows the result 



that electrons and holes are excitonically bound. In our 
trapping configuration it is only possible to confine IX, 
free electrons are repelled. For free charge carriers or an 
electron-hole plasma different dependence of the PL line 
is reported, reflecting quantized Landau levels [5J [TU] . 

The thermal de Broglie wavelength AdB of a bosonic 
particle in two dimensions with a quadratic dispersion 



relation is given by AdB 



where m is the total 



V2rmrkT 

mass of the exciton (sum of electron and hole mass) [11] . 
The interexcitonic distance dix in a two-dimensional ex- 
citon lattice hexagonal arrangement is given by dix = 

In the correlated many-body state the thermal 
de Broglie wavelength (« 75 nm at 2 K) is larger than the 
interexcitonic distance drv 57 nm at n = 4 x 10^° 
and dix is similar to twice the effective Bohr radius ~ 
20 nm. 



DIAMAGNETIC SHIFT OF THE EXCITON 
ENERGY IN MAGNETIC FIELDS AND 
RELEVANT LENGTH SCALES 



ANALYSIS OF THE RED TAIL OF THE 
SPECTRA 



A magnetic field B applied in Faraday configuration 
perpendicular to the QW plane causes a diamagnetic 

quadratical shift of the electrostatically trapped IX en- 

2 2 

ergy by A£^dia — g^^^ where is the effective Bohr 
radius, ^ the reduced mass of the exciton [71 [H]. In 



Fig. S5 the quadratic behavior of the PL maximum en- 



ergy, detected from the center of the trap, is shown as a 
function of the magnetic field. Fitting the exciton dia- 
magnetic shift we find a corrected effective Bohr radius 
of r* = 20 nm [S]. If the hole is strongly localized, it 
is necessary to use the electron mass instead of the re- 
duced mass fi which results in an effective Bohr radius 
of r** = 16 nm. The quadratic diamagnetic shift shows 



1.398 




2 4 

Magnetic field (T) 



FIG. S5: Magnetic field dependence. The PL maximum 
detected in the center of the trap is plotted as a function of 
the magnetic field applied perpendicular to the QW plane. 
In addition, a quadratic fit for the diamagnetic shift is shown 
from which we find an effective Bohr radius rc — 20 nm [Vg 
= 0.70 V, Vt = 0.32 V, and Vs = 0.35 V, TLatticc = 256 mK, 
-BLaser ~ 1.4289 eV, Pbaser = 3.5 /uW circular polarization]. 



The analysis of the red tail of the photoluminescence 
spectra clearly showed the occurrence of different power 
laws with distinct exponents. The first exponent with a 
value of about -0.1 can be related to a low-energy exci- 
tation of up to 1.5 meV and reflects slow temporal be- 
havior, possibly induced by thermal broadening. The 
second exponent shows values of about -0.8 describing 
a medium energy excitation of the IX ensemble in the 
energy interval between 1.5 and 4meV. This power law 
rather abruptly disappears at TLattice ^ 5K [Fig. 2(e) 
in the main text] and we relate it to a correlated many- 
body phase of a liquid of parallel oriented dipolar exci- 
tons as theoretically predicted in Ref . [5] . This correlated 
phase appears only at low temperatures and high enough 
trapped IX densities. The third exponent of about -3 is 
related to high-energy {Eq — E > 4 meV) excitations of 
the IX ensemble on a fast time scale. Most likely, phonon 
emission on a picosecond time scale and shake-up pro- 
cesses, in which the recombining exciton transfers part 
of its energy to a collective excitation of its interacting 
neighboring excitons, are responsible for this behavior. 

In order to determine the exponent of the different 
power laws I{E) ~ (Eq — i?)~l"l in the red tails of the 



PL lineshape in the log- log plots [Fig. S6[a), (b)] we em- 

In Fig. [S6] this 
In the 



ployed the following numerical method, 
method is shown exemplarily for one spectrum, 
first step the first derivative of the logarithmized data 
is calculated. Simply taking the difference between two 
neighboring data points would yield a very noisy result. 
For a fixed energy Eq — E a, linear fit is done over ±i 
neighboring data points. The slope obtained by the lin- 
ear fit over the 2i-f 1 data points gives the value of the first 
derivative for this energy. This is repeated successively 
for each energy point. The result is shown in Fig. S6 [c). 
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The second derivative is calculated in the same way from 
data points of the first derivative. All data pairs whose 
values of the second derivative are in a range of ±0.5 are 
considered as having a constant slope. These data points 



are between the two horizontal lines in Fig. S6 'd) . If the 
slopes of two or more data points differ less than 0.1 for 
the first slope and 0.5 for all other slopes these slopes and 
their corresponding energies are averaged. The previous 
values provide the best results but small variations do 
not change the exponents of the power laws. In contrast 



to Fig. 2(d) in the main text in Fig. S6 a) the logarith- 



mized data are plotted to illustrate that in our numerical 
method a linear fit of the logarithmized data is done. 

We used a second method as a consistency check by 
performing a linear fit over ±i neighboring data points. 
If none of the 2i + 1 residues is larger than a predefined 
value, these data points are considered as having a con- 
stant slope. This procedure is also repeated successively 
for each energy point, and data points with similar slopes 
are averaged in the same way as described above. Both 
methods yield comparable results. However, the first 
method turned out to be more general since the signal- 
to-noise ratios of our spectra differ substantially and the 
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FIG. S6: Analysis of the red tail of the spectra in a log- 
log diagram, (a) log-log plot of one typical spectrum. The 
slopes are determined by the numerical method described in 
the text [Vg = 0.70 V, Vt = 0.32 V, and Vs = 0.34 V, PLasor 
= 1.69 /^W, TLattice = 248 mK, i = 5]. (b) Original spectrum 
with power laws determined by the slopes, (c) First deriva- 
tive of the logarithmized data, (d) Second derivative of the 
logarithmized data. The values between the two horizontal 
lines represent data points which were used to calculate the 
slopes. 



residues depend strongly on the signal-to-noise ratio. 



INTERFERENCE SUBSTRUCTURE 
SUPERIMPOSED ON THE PL LINESHAPE 



Photoluminescence spectra having many photon 
counts and correspondingly large signal-to-noise ratio 
show a pronounced substructure in the form of energeti- 



cally equidistant shoulders. In Fig. S6 ^b) this shoulders 
are marked with arrows. While the spectra shift with 
changing IX density the substructure is fixed on an abso- 
lute energy scale. We identify the energetic periodicity of 
the substructure as an interference effect inside the het- 
erostructure, where a cavity is formed between the top 
and the bottom surface of the sample. A similar interfer- 
ence effect has been reported in [H]. We found that the 
slopes are influenced by the substructure and the slope 
starts to oscillate. Nevertheless, the data at low densi- 
ties yield a convincing evidence that the power law with 
-|a| « -0.8 is an intrinsic and characteristic feature and is 
not caused by such interference effects. All temperature- 
dependent measurements show that the slope slightly 
lower than 1 disappears for temperatures T ^ 5 K. This 
suggests a breakdown of the correlation effect and cannot 
be explained with an interference effect. In log-log-plots 
of spectra with normal signal-to-noise ratios we observe 
no pinning of the boundary point positions on the inter- 
ference substructure. We explain this with the absence 
of the substructure due to the lower signal-to-noise ratio. 



INDIRECT EXCITON GAS TEMPERATURE 



We use two different methods to determine the appar- 
ent exciton gas temperature Tapp [see also Fig. 2(c) and 
Fig. 3(c) in the main text] from the blue tail of the PL 
spectrum. In the first method the exciton gas tempera- 
ture is extracted from the exponential Boltzmann-like be- 
havior of the blue tail. Therefore the natural logarithm of 
the PL intensity is calculated and the slope of the linear 
regime is determined with the same numerical method as 
described above. One slope obtained with this method 



is shown in Fig. S7 [a,) together with the logarithmically 



plotted spectrum. Using logarithmic data for the fit has 
the advantage that data with low intensity away from 
the PL maximum are weighted more strongly compared 
to an exponential regression of the original data. In this 
regime the assumption of a Maxwell-Boltzmann distribu- 
tion is approximately valid. The exciton gas temperature 
can be calculated directly from this slope. Our second 
approach is to multiply a Gaussian distribution with a 
Fermi-Dirac distribution. The Gaussian distribution is 
introduced to reflect broadening of the PL line by (i) the 
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energy dependent density of states of the interacting IX 
dipoles in the trap (ii) the random disorder potential in 
the CDQW and (iii) effects of time-dependent spectral 
diffusion as the displayed spectra are accumulated over 
approximately 1 minute. The Fermi-Dirac distribution is 
used as a cutoff function, reflecting the thermal smearing 
of the filling of the trap with the dipolar exciton liq- 
uid. In combination the Gaussian distribution describes 
the red tail of the spectrum while the Fermi-Dirac dis- 
tribution describes the temperature-broadened steep de- 
crease on the blue tail of the spectrum. In the case of the 
bosonic dipolar liquid the strong spatial dipolar repulsion 
replaces the role of the Pauli principle in the Fermi-Dirac 
distribution. One spectrum together with its fit is shown 
in Fig. S7'b). Although the assumption of a Fermi sea 



Boltzmann fit. Both fits show a constant offset to the 



can only be considered approximatively since we have an 
IX ensemble with hybridized electrons and no free elec- 
trons, the results of this method are in good agreement 
with the exciton gas temperatures obtained by the first 
method. The temperatures extracted by the Fermi-Dirac 
distribution scatter less than these using the Maxwell- 



lattice temperature [green line in Fig. S7'c)] which at 
least in part can result from the broadening contributions 
discussed above which persist to lower lattice tempera- 
tures. Additional temperature independent broadening 
of the blue tail of the excitonic emission, not included in 
the modeling, is likely to result from scattering events or 
shake-up processes and may be an essential contribution 
to this constant offset of the exciton gas temperature ex- 
tracted from the fits. Hence Tapp should be considered 
as an upper limit of the IX temperature whereas the real 
IX temperature is likely to be considerably closer to the 
lattice temperature. For a lattice temperature TLatticc ^ 
1 to 2 K, we observe a saturation of the Tapp [Fig. S7 c)] 



corresponding to an upper bound of the apparent exci- 
ton gas temperature at Tapp = 3 K. Fig. S7 [d) shows the 



saturation temperatures of several temperature sweeps 
with different excitation powers. This plot indicates that 
the saturation temperature decreases exponentially with 
decreasing excitation power. 




^Lattice ^Laser t^*^) 



FIG. S7: Extraction of an apparent exciton gas tem- 
perature, (a) Natural logarithm plot of the blue side of one 
spectrum. The red line is obtained by the numerical method 
described in the text [Vb = 0.70 V, Vr = 0.32 V, and Vs = 

0.35V, PLascr = lAflW, TLatticc = 285 mK, i = 5]. (b) PL 

spectrum from the same data as in (a). The red curve repre- 
sents a "Gauss-Fermi-Fit", which is described in the text, (c) 
Fit results of a temperature dependent measurement. Both 
methods are in good agreement, (d) Saturation temperatures 
of Tapp for different excitation powers extracted from different 
temperature sweeps like those displayed in (c) [Vg ~ 0.70 V, 
Vt = 0.32 V, and Vs = 0.35 V]. 
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